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Abstract 

The femtosecond magnetization dynamics of a thin cobalt film excited with ultrashort laser 
pulses has been studied using two complementary pump-probe techniques, namely spin-, energy- 
and time-resolved photoemission and time-resolved magneto-optical Kerr effect. Combining the 
two methods it is possible to identify the microscopic electron spin-flip mechanisms responsible 
for the ultrafast macroscopic magnetization dynamics of the cobalt film. In particular, we show 
that electron-magnon excitation does not affect the overall magnetization even though it is an 
efficient spin-flip channel on the sub-200fs timescale. Instead we find experimental evidence for 
the relevance of Elliott- Yafet type spin-flip processes for the ultrafast demagnetization taking place 
on a time scale of 300 fs. 

PACS numbers: 75.70.Ak, 78.47.+p, 79.60.Dp 
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Understanding the microscopic processes governing magnetization dynamics on the fem- 
tosecond time scale is one of the most challenging and interesting issues in today's condensed 
matter physics. There is nowadays much experimental evidence that the excitation of a fer- 
romagnetic film with an ultrashort laserpulse gives rise to a demagnetization on a time 




scale of a few hundred femtoseconds jjj, |2|, |3, y, |5|, |6|] . However, the microscopic processes 
responsible for these ultrafast changes in the macroscopic magnetism of the ferromagnet 
have not been unambiguously identified yet. A complete understanding of such processes 
would be extremely helpful for technical applications like ultrafast magnetization-switching 
schemes for read-write processes. Most of the recent experimental approaches applied in 
this field employed time-resolved magneto-optical strategies Such experiments 

access the ensemble spin dynamics taking place on the femtosecond time scale. On the other 
hand, spin- and time-resolved photoemission addresses the dynamics of single electron spin- 
flip processes, which are responsible - from a microscopic point of view - for the observed 
changes in the sample magnetization. Bringing together these two experimental methods 
has, therefore, the potential for an efficient access to the microscopic view of the ultrafast 
dynamics of macroscopic magnetism. 

In this letter such a complementary approach is presented for the first time and is used 
to obtain a more comprehensive understanding of the ultrafast magnetization dynamics of 
a thin cobalt film. Time-resolved magneto-optical Kerr effect (TR-MOKE), as well as spin-, 
energy- and time-resolved two-photon photoemission (SETR-2PPE) are used to address re- 
spectively the macroscopic and microscopic dynamical response of the ferromagnet induced 
by ultrashort optical excitation. SETR-2PPE probes the time evolution of the spin polariza- 
tion of single optically excited electrons at a time-resolution down to a few femtoseconds. In 
comparison with TR-MOKE data we find that characteristic features in the observed com- 
plex time dependence of the SETR-2PPE signal can be directly correlated to changes in the 
magnetization of the sample. In this way, we are able to draw conclusions about the micro- 
scopic spin-flip mechanisms leading to the observed macroscopic ultrafast demagnetization 
of the thin Co film. 

There are three main mechanisms through which an excited electron can flip its spin 
in a ferromagnetic metal: (i) a Stoner excitation (ii) an inelastic electron-spinwave scat- 
tering event (magnon excitation), or (iii) an Elliott- Yafet type of spin-flip scattering with 
impurities or phonons [3, Whereas for excitation energies above the vacuum energy 



2 



(E — Ep > $, $=work function), channels (i) and (ii) have been identified within static 

nn 

SPEELS measurements |9|,110J], there exists no direct experimental access to optically ex- 
cited electrons which are most relevant for laser induced ultrafast demagnetization. Thus, 
the excitation regime essential for ultrafast laser demagnetization (E — Ep < $) remains 
to be explored experimentally. Theoretically, short wavelength spinwave excitations have 



been studied in 
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12j. According to these studies, spinwave scattering dominates Stoner 



excitation processes in the very low energy regime (roughly (E — Ep) < 3eV). Moreover, 
first-principles GW and T-matrix calculations including spin-flip scattering 13] have shown 
that the coupling to spinwaves is an important inelastic decay channel of excited electrons 
in itinerant ferromagnets, especially for energies (E — Ep) < 1.5 eV. On the other hand, a 
detailed analysis of TR-MOKE experiments has provided some evidence that Elliott- Yafet 
processes have to be considered for a correct interpretation of femtosecond magnetization 
data |l4|. In the following we will show that electron- magnon interaction indeed contributes 
significantly to the spin dynamics of single electrons excited close to the Fermi level (Ep) 
in a Co thin film. However we also find that this process cannot account for the observed 
ultrafast demagnetization. Instead we identify clear support for the Elliott- Yafet scenario. 

The detailed experimental setup for SETR-2PPE has been described elsewhere [l^- 
Briefly, we use the second harmonic of a femtosecond Ti: sapphire laser with photon en- 
ergy hv = 3.1 eV. The beam is divided into two equally intense, collinear and orthogonally 
polarized beams with an adjustable time delay. The pump photons (here, s-polarized) excite 
electrons from their ground state into an unoccupied intermediate state. The delayed probe 
photons (p-polarized) excite the electrons above the vacuum level. The angle of incidence 
of the laser on the sample is 45°; the fluence of the pump and probe pulses is 3/^J/cm 2 . The 
photoemitted electrons are analyzed in normal emission geometry in an ultrahigh vacuum 
system equipped with a commercial cylindrical sector energy analyzer (Focus CSA 300) and 
a spin detector (Focus SPLEED). 

Co films with a thickness of 10 nm are grown quasi-epitaxially on a Cu(001) substrate, 
and possess in-plane magnetization easy axis along the (110) direction of the Cu crystal. 
For the SETR-2PPE experiments the Co films are remanently magnetized by an external 
field applied along the easy axis direction. The spin detector measures the projection of the 
spin polarization along the easy axis: P = (N + — N~)/(N + + N~), where N + (N~) is the 
number of detected majority (minority) electrons. 
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After the SETR-2PPE experiments, the same Co film has been investigated ex-situ by 
TR-MOKE. A two-color pump-probe experiment has been used, where the pump pulses are 
generated by a multipass amplifier system operating at a central wavelength of 800 nm and 
delivering 50 fs pulses with 1.6mJ/pulse at a repetition rate of 1kHz. The probe pulses are 
obtained by frequency doubling a small fraction 1/10) of the pump beam in a /3-barium- 
borate crystal. The bichromatic approach enables an efficient suppression of bleaching effects 
which may otherwise lead to a misinterpretation of the MOKE data j^. The experiments 
were performed by recording, for every pump-probe delay, full hysteresis loops in longitudinal 
configuration. In this way the shape of the hysteresis can be used as a monitor of the sample 
quality during measurements. The external field vector is applied parallel to the in-plane 
easy axis direction. The inset of Figure |3] shows two exemplary hysteresis loops recorded 
respectively without the pump pulse (black curve) and 675 fs after it (grey curve). 

In a time- and energy-resolved 2PPE measurement the recorded quantity (cross- 
correlation curve) is the number of excited electrons at a certain energy as a function of 
the time delay between the pump and the probe pulse (population decay). Figure [T] (a) and 
(b) show the cross-correlation curve of an intermediate state with an energy of 2.8 eV and 
0.4 eV above Ep, respectively. In case the 2PPE signal is governed exclusively by an expo- 
nential population decay of the intermediate state, the experimental data can be fitted by 
a rate-equation model, where the lifetime r of the selected intermediate level is the only fit 
parameter. Corresponding calculations to simulate our experimental data were performed 
using the formalism described in in the limit of rapid dephasing between the electronic 
excitation and the driving laser field. The dashed lines in Figure d correspond to such a fit. 
Whereas the data recorded at 2.8 eV intermediate energy (a) are almost perfectly described 
by this model, we find for the lower energy (b) rather significant deviations at temporal 
delays between 50 and 200 fs. In order to quantify these deviations, we evaluate the refilling 
R = (y m — Yf)/(Y m + Yf), where Y m is the measured 2PPE yield, and Yf the result of 
the rate-equation model. The obtained curves are shown on the right scale (open circles). 
Clearly, at 2.8 eV intermediate energy R is almost zero, while at 0.4 eV it is significantly 
different from zero up to 200 fs delay, and reaches its maximum value of ~ 1 % at 120 fs. The 
explanation for this additional contribution to the 2PPE signal is that secondary electrons 
scatter into the probed intermediate level and lead to a repopulation (refilling) of the state 
itself (see Figure 12(c)). From the values of R we can conclude that in our case refilling elec- 
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trons constitute at maximum 1 % of the electron population at the considered intermediate 
level. There are different relaxation mechanisms resulting in refilling electrons; the most 
relevant is the scattering of an optically excited (hot) electron in a higher excited state with 



a cold electron from below Ep. Further details on refilling can be found in 

Let us now concentrate on the time evolution of the spin polarization of the refilling 
electrons in the intermediate state located 0.4 eV above Ep as detected by SETR-2PPE. 
Our aim is to identify which kind of spin-flip scattering events take place during the decay 
process as schematically shown in Figure ^(c). Figure El shows the dependence of the in- 
plane polarization of electrons in the intermediate state at this energy. In particular, we 
plot the normalized polarization changes, defined as AP = (P — -P m in)/(-Po — -Pmm), where P 
is the polarization at zero delay and P m i n the minimum of the polarization over the studied 
delay range. AP exhibits a complex time-dependence with two maxima at 30 fs and 120 fs, 
and two minima at 75 fs and 350 fs. 

To separate different contributions to the spin- and energy-resolved electron dynamics, 
we use a numerical model which phenomenologically includes carrier Boltzmann scattering 
together with carrier transitions due to pump and probe pulses. This model shows which 
features of the experimental results can be explained without the existence of magnetic cor- 
relations (such as magnon excitation). More precisely, we solve dynamical, energy-resolved 
equations for the electronic occupation n s (E) of single-particle states with spin s and energy 
E of the form 



^n s {E) = - TT{E)n s {E) + T™{E) [1 - n s {E)\ 
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where the last term on the right-hand side describes electronic transitions due to the pump 
and probe laser fields and includes the density of states which is taken from experiment. 
The in and out-scattering rates in Eq. (£[]) for weak nonequilibrium situations are determined 
from the measured spin-dependent electron lifetimes using T° nt (E) = —b(—E)/r s {E), where 
b(E) = [exp(PE) — l]^ 1 is the Bose function, together with the detailed-balance relation 
r° ut (£') = exp(l3E)r™(E). Using these relations implies that the dynamics obtained from 
Eq. (^) describes transitions involving scattering events with both minority and majority 
electrons, but does not account for the full correlated electron dynamics. In particular, 



our approach can describe the influence of the dynamical single-particle occupation on the 
experiment, which includes effects such as cascade electronic scattering events. 

The result of the simulation for delays up to 200 fs is shown in the inset of Figure|21 (dashed 
line), together with the measured data. Clearly, a good agreement between experiment and 
calculations is achieved only for the first 100 fs. Here, the spin polarization first increases 
from its value at zero delay to a maximum close to 30 fs. This is due to the longer lifetime of 
majority electrons compared to the lifetime of minority electrons [17]. This difference implies 
that minority electrons depopulate the intermediate state at a faster rate than majority 
electrons, leading to a transient increase in the spin polarization (spin-filter effect). Between 
30 and 75 fs the spin-filter effect starts to become irrelevant and the polarization decreases 
to the background value given by the spin-polarization excited by the pump and probe 
pulse separately. The successive increase of the spin polarization in the experiment, with a 
maximum around 120 fs, is not reproduced by the simulation. At this time delay we have a 
maximum of the refilling curve R in Figured fb). In order to interpret this behavior correctly, 
let us first observe that the calculations in |lj| predict that minority electrons exhibit a 
higher cross-section for magnon excitation compared to majority electrons. Combining this 
results with the fact that our numerical simulation does not include magnon excitation, we 
can conclude that this process causes the measured increase of the spin polarization around 
120 fs. In other words, the intermediate state at 0.4 eV above Ep is refilled by secondary 
electrons originating from energetically higher-lying minority electrons (with energy up to 
3.1 eV above Ep): those electrons have a higher probability than majority electrons to scatter 
inelastically with a magnon and flip thereby their spin, turning into majority electrons and 
thus causing an increase in the measured spin polarization (see Figured ( c ))- 

We now discuss the behavior of the spin polarization for time delays between approxi- 
mately 150 and 1000 fs. It shows a drop with a minimum around 350 fs and a successive 
recovery. For such time delays, the number of pump laser excited electrons is so small that 
the observed behavior must be due to a change in the magnetic state of the electrons within 
the Fermi see. Indeed, after optical excitation the excited electrons and the produced sec- 
ondary electrons thermalize to a hot Fermi-Dirac distribution on such a time scale. This 
implies that the feature observed in Figure El is related to a process in which the total mag- 
netization is changing, and should be observed also with TR-MOKE. Figure HU shows the 
TR-MOKE signal recorded for 10mJ/cm 2 (triangles) and 2mJ/cm 2 (squares) pump fluence. 
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The curves are obtained by plotting, for a given pump-probe delay, the Kerr amplitude of 
the corresponding hysteresis loop, normalized by the value at zero delay. They show an 
ultrafast drop followed by a recovery. By reducing the pump fluence the minimum of the 
curve shifts to lower time delays. For the lower fluence the minimum is reached at approx- 
imately 400 fs, close to the minimum observed in the SETR-2PPE experiment. This fact 
gives a clear indication that the mechanism responsible for the reduction of the polarization 
in the SETR-2PPE experiment is the same one causing a reduction of the TR-MOKE signal, 
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i.e. it is a process in which the total magnetization is changing. Koopmans et al. 
have argued that the ultrafast demagnetization observed in TR-MOKE experiments can be 
microscopically explained by assuming Elliott- Yafet type of spin-flip scattering of electrons 
with impurities or phonons. According to the same authors, such process leads to a decrease 
of the minimum of the TR-MOKE curves with the pump fluence. Such a behavior was in- 

n 

deed observed in [19||, though, at lower fluencies and thereby with more subtle variations. 
Our TR-MOKE results give a valuable support to such interpretation and moreover, through 
the comparison with SETR-2PPE provide the experimental proof of the dominating role of 
Elliott- Yafet scattering in a microscopic description of ultrafast demagnetization. 

In conclusion we have successfully combined two experimental techniques, one (TR- 
MOKE) giving an insight into the ensemble femtosecond spin dynamics after laser excitation, 
the other (SETR-2PPE) providing additional information on the microscopic spin-flip pro- 
cesses responsible for the observed magnetization dynamics. We were able to conclude that 
even if microscopically spinwave excitation plays an important role in the ultrafast spin 
dependent electron dynamics, the ultrafast demagnetization of the studied thin Co film is 
mainly due to Elliott- Yafet type spin-flip processes. 
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FIG. 1: M. Cinchetti, Physical Review Letters 
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FIG. 2: M. Cinchetti, Physical Review Letters 
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CAPTIONS 



FIG. 1: (a) and (b). Left scale: 2PPE yield from an intermediate state with energy of 
2.8 eV (a) and 0.4 eV (b) as a function of the delay time between the pump and the 
probe pulse. The fit of the experimental curves with a rate-equation model is plotted 
with dashed lines. Right scale: Refilling parameter R as derived from the measured 
curves and the fitting curves. 

(c) Schematic picture of a possible spin-dependent refilling process of the intermediate 
state at 0.4 eV, as discussed in the text. 

FIG. 2: Time dependence of the spin polarization of electrons in the intermediate state 
at 0.4 eV above Ep as detected by SETR-2PPE. Inset: a zoomed picture of the same 
curve for time delays between Ofs and 200 fs and the result of the numerical model 
discussed in the text (dashed line). The continuous line (obtained by smoothing the 
original data) is depicted as a guide to the eye for the reader. 

FIG. 3: Results of TR-MOKE measurements on the same Co sample of Figure |2j taken 
with a pump fluence of 10mJ/cm 2 (triangles) and 2mJ/cm 2 (squares). Inset: two 
exemplary hysteresis loops recorded respectively without the pump pulse (black curve) 
and 675 fs after it (grey curve). The pump power was 10mJ/cm 2 . 
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